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@ Object-based irregular-grid volume rendering. 

@ Method and apparatus for rendering data for display. The data represents a model space and the method 
includes the following steps. A first step associates individual points In the model space with individual vertices 
(V) of an irregular grid cell enclosed by faces. Each of the points has an associated luminosity value and an 
associated opacity value. A next step orders the faces of the grid cells. The faces are preferably ordered by 
depth referenced to a view plane. A further step processes individual faces for determining, for all viewing rays 
that project from pixels on the view plane and that intersect the face being processed, an amount of light aniving 
at the pixel. The amount of light arriving at the pixel is a function of both the luminosity and the opacity of a grid 
cell that lies adjacent to the face and through which the viewing ray passes. 
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This invention relates generally to Image display apparatus and method and. in particular, to an object- 
based method and apparatus for rendering volumetric data that is represented within an irregular grid. 

Volume rendering Is a known technique for visualizing three-dirnensional space-filling, or volumetric, 
data. In general, volume rendering assigns to each point within the volume of data an opacity and a 
5 luminosity. The opacity and luminosity are assigned as some function of the data value at that point. 
Subsequently, there is produced an image of the resulting imaginary solid translucent object. By proper 
choices of colors, opacities, and viewing directions a visualization and, hence, an understanding of the data 
may be obtained. 

Traditionally, most applications have operated with data sampled onto a uniform or regular grid, such as 
10 medical computerized tomography Images. Volume rendering of regular grids is readily accomplished by 
compositing successive layers of volume elements, or voxels, to obtain an image. However, many scientific 
simulations are accomplished with Irregular grids, in which the volume of interest is divided into a space- 
filling mesh of cells where the data values are computed only at the vertices where the cells meet. Data 
values within a cell are interpolated from the data values at the vertices of the cetL As one example, and as 
76 illustrated in Fig. 9, an irregular grid comprised of cells 1 may be employed for modelling airflow over a 
curved surface of an aircraft wing 2. The grid is made irregular by conforming the grid to the surface under 
study. Of course, this is but one example of an irregular grid structure. 

The relatively simple compositing method does not apply to rendering such irregular grid data because 
the thickness of each cell, as viewed from each pixel at the image plane, must be taken into account in 
20 forming the image. Rendering such irregular-grid volume data currently requires either resampling the data 
on a regular grid or using a ray tracing technique. 

Unfortunately, neither of these conventional approaches provides an optimum rendering solution. 
Resampling and rendering the data onto a regular grid typically greatly increases the number of data points 
to be processed. Ray tracing is also a computationally expensive technique in which imaginary viewing rays 
25 are projected out from the image plane. The luminosity and opacity are computed at each cell, taking cell 
thickness Into account, and are then integrated along each ray. 

The following U.S. patents are cited for showing conventional regular grid volume rendering systems. 
U.S. Patent 4,866,612, issued September 12, 1989, entitled "Imaging System Including Means to 
Provide A Display Image with Shadows and Three Dimensional Perspectives by a Combination of Voxel and 
30 Surface Method" to Takagi et al. relates to displaying surfaces derived from a volumetric regular grid. 

U.S. Patent 4,827.413, Issued May 2, 1989, entitled "Modified Back-To-Front Three Dimensional 
Reconstruction Algorithm" to Baldwin et al. describes a method for displaying opaque solids defined by a 
volumetric regular grid. 

U.S. Patent 4.719.585. Issued January 12, 1988, entitled "Dividing Cubes System and Method for the 
35 Display of Surface Structures Contained within the Interior Region of a Solid Body" to Cllne et al. renders a 
surface defined by a threshold value on a regular volumetric (parellelopiped) grid. 

U.S. Patent 4,835,688, Issued May 30, 1989, entitled "Three-Dlmensional Image Processing Apparatus" 
to Kimura renders a surface defined by a threshold value on a regular volumetric grid. 

U.S. Patent 4.821.210 issued April 11, 1989. entitled "Fast Display of Three-Dimensional Images" to 
40 Rumbaugh extracts a surface In the form of a set of triangles directly from volume elements. Specific 
mention is made of extracting polygonal surfaces froni cubes. 

U.S. Patent 4.729,098, Issued March 1, 1988, entitled "System and Method Employing Nonlinear 
Interpolation for the Display of Surface Structures Contained Within the Interior Region of a Solid Body" to 
Cllne et al. extracts a surface in the form of a set of triangles directly from the volume elements. This patent 
45 also considers only regular grids and has disclosure directed to extracting polygonal surfaces. 

An article entitled "Chem-Ray: A Molecular Graphics Program Featuring an Umbra and Penumbra 
shadowing Routine" by J. W. Lauher, J. Mol. Graphics, 1990, Vol. 8, March pp. 34-38 uses a ray tracing 
technique to render shadowed molecular graphics models. 

Also of interest in this area is U.S. Patent 4.835,712. issued May 30. 1989, entitled "Methods and 
50 Apparatus for Imaging Volume Data with Shading" to Drebin et al.. This patent describes a method of 
assigning material opacities to volume elements so as to make boundaries within volumetric data visible. 
Drebln's teaching makes mention of storing data In arrays of voxels. Implying regular grids. Drebin et al. 
describe a method for computing opacities and luminosities at each point In the volumetric data, but do not 
specify how the resulting opacities and luminosities are to be rendered for display. 
55 It is an object of the invention to provide an efficient method and apparatus for rendering volumetric 

data. 

The invention especially provides method and apparatus for rendering volumetric data defined at 
sample points on an irregular grid, such as a mesh of tetrahedra, deformed rectangles, or other non- 
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rectangular volume elements, and In particular an object-based, as opposed to a ray-based, method and 
apparatus for rendering volumetric data defined at sample points on an irregular grid. 

Preferably the invention further provides an efficient method and apparatus for rendering volumetric 
data as if the data represented a solid in space having a specific opacity and luminosity at each point, as 
5 opposed to only rendering surfaces derived from volumetric data or to rendering the volumetric data as if it 
were an opaque solid. 

This and other objects are advantageously solved basically by applying the features laid down in the 
independent claims. 

The object of the invention is realized by an object-based method and apparatus for volume rendering 

10 three-dimensional data defined on an irregular grid, such as scientific or engineering simulation results. In 
accordance with the invention there is provided a rendering technique in which each face, or border 
between grid cells, is considered in turn. The visual effect of a cell that lies in front of or in back of the 
considered face is added to the image. This technique provides a significant speed advantage over the 
conventional techniques of resampling and ray tracing. 

75 In accordance with a method of the invention, and apparatus for accomplishing the method, there is 
described a method of rendering data, representing a model space, for display. The method includes the 
following steps. A first step associates individual points in the model space with individual vertices of an 
irregular grid cell enclosed by faces. Each of the points has an associated luminosity value and an 
associated opacity value. A next step orders the faces of the grid cells. The faces are preferably ordered by 

20 depth referenced to a view plane. A further step processes individual faces for determining, for all viewing 
rays that project from pixels on the view plane and that intersect the face being processed, an amount of 
light arriving at the pixel. The amount of light arriving at the pixel is a function of both the luminosity and the 
opacity of a grid cell that lies adjacent to the face and through which the viewing ray passes. 

The step of associating includes a step of determining an x, y and z coordinate of the point relative to 

25 the view plane and also includes a step of expressing the luminosity in terms of a plurality of colors. 

The step of processing may be accomplished by processing the ordered faces in order of increasing 
depth from the view plane or in order of decreasing depth from the view plane. 



The above set forth and other features of the invention are made more apparent in the ensuing detailed 
description of the invention when read in conjunction with the attached drawing, wherein: 



30 


Fig. 1 


is a simplified block diagram of a graphics rendering system constructed and operated 






in accordance with the invention: 




Fig. 2 


shows in greater detail the relationship of a point list and a vertex list that are 






maintained in the memory of the system of Fig, 1; 




Fig. 3 


shows in greater detail the organization of a face list that is also maintained in the 


35 




memory of the system of Fig. 1 ; 




Fig. 4 


depicts a viewing ray, a viewing plane, and their relationship to volumetric data to be 






rendered; 




Fig. 5 


depicts a face of a tetrahedral grid cell having vertices VI . V2, and V3; 




Figs. 6 and 7 


depict a ray tracing technique of the prior art; 


40 


Fig. 8 


depicts a viewing ray (r) passing through a plurality of irregular, non-rectangular, grid 






cells; and 




Fig. 9 


is a diagram that illustrates an irregular grid, of the type processed by the invention. 






that is applied to a surface of a wing for modelling an airflow pattern. 



Referring first to Fig. 1 there is depicted a rendering system 10 that is constructed and operated in 
45 accordance with the invention. System 10 includes a digital bus 12 that couples together a data processor 
14. a memory 16, a frame buffer 18 that is coupled to a display means 20, and a display controller 22 
having a user entry device, such as a keyboard 24. Preferably the display 20 is a high resolution color 
graphics display. The processor 14 processes volumetric data that is proved thereto to create a displayable 
image in an image buffer (IB). The contents of IB are moved across the bus 12 to the frame buffer 18 for 
50 display. In accordance with the invention the volumetric data provided to the processor 14 is referenced to 
an array of irregular, non-rectangular grid cells. By example, the grid cells may be tetrahedral, rectangular, 
or of any three dimensional polyhedral shape selected for representing the underlying volumetric data to be 
rendered and visualized. To this end the* processor 14 employs a first Input list (LI) of points or vertices of 
the grid cells, a second Input list (L2) of elements that specify the vertices of the grid cells, and a third 
55 internally generated list (L3) of ordered faces of the irregular grid cells. 

As employed herein it is assumed that a solid to be rendered has defined at each point a luminosity, 
expressed as quantity of light emitted per unit of material, and an opacity (or optical density), expressed as 
a proportion of light absorbed per distance travelled. The opacity is generally an Intrinsic property of the 

3 
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material, related In some way to the volumetric data. The luminosity may be intrinsic, that is derived from 
the volumetric data; or extrinsic, that is due to external lighting or internal shadows or scattering. 

Before continuing a further discussion of the system shown in Rg» 1 a description will first be provided 
of the underlying mathematical relationships implemented by the system in rendering volumetric data. 

5 Referring briefly to Fig. 4 an amount of light arriving at any point in a viewing plane is determined by 

integrating the contributions from each point along a viewing ray between a near point zn and a far point zp. 
If L(2) and D(2) represent the luminosity and opacity, respectively, at a point z along the ray; the light 
emitted from each point z on the ray toward the viewing plane is attenuated by passage through the 
material from the point z on the ray to the point zn. The amount of this attenuation is given by the 

70 transparency T: 



75 
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T(zN,z)*^i^exp|- D{C)d; 

'zN 



\Jz 



(1) 



It should be noted that while D(2) indicates the proportion of light absorbed per distance travelled at the 
point 2; T(2n.2) indicates the proportion of light transmitted in travelling from z to Zn. To a first approxima- 
20 tion, 

T(2,2 + A2)«1-D(z)A2. (2) 

The total amount of light arriving at the viewing plane is given by integrating the light L(z)dz emitted at each 
25 point 2, attenuated by the transparency T(2n. z) from z to 2n. in general, the amount of light emitted toward 
the viewing plane by a segment of a viewing ray from Zu to 2f, or the brightness of that segment, is given 
by the integral 



rzF 

B(zN, 2f) 4?^ L (z) T(ZN, z) dz 





(3) 



T and B can be shown to exhibit simple but useful properties. A first property is that the total 
40 transparency of two adjacent segments is the product of their transparencies: 

T(a,c) = T(a,b)T(b.c), where a ^ b ^ c. (4) 

A second useful property is that the brightness of a segment may be computed from the brightness of 
45 two sub-segments according to: 

B(a,c) = B(a.b) + T(a.b)B(a.c). (5) 

where a ^ b S c. 

50 By applying the properties of T and B from equations 4 and 5 the brightness (B) can be rewritten in 
terms of the transparencies and brightnesses of each of the segments (Zj. Zj+i) as: 

Computing B is accomplished by: enumerating the intersections between the viewing rays, originating 
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at each viewing plane pixel, and the faces; at each intersection interpolating to obtain the lunninosity and 
opacity; for each segnnent. approximating Bj and Tj; and accumulating the light for each viewing ray 
according to the sum in equation (6). 

In accordance with an aspect of the invention this calculation of B is carried out most efficiently by an 
object-based technique thai enumerates faces in front-to-back or back-to-front order, and for each face 
enumerating that face's intersection with viewing rays. 

Let X denote the accumulated light for any given pixel. For front-to-back order. It is also necessary to 
accumulate t, the total transparency of the segments processed thus far. These variables are initialized as: 

X — 0, and 
T — 1 . 

It should be noted that in the case of a color display there is provided a separate X, and a separate r if 
the material absorbs colors selectively, for each of primary, display colors red, green, and blue (RGB). 

For each segment (2|, Zj+i) there is computed an estimate Bi for Bj and an estimate tt for t( in 
accordance with 

AZj "^21 + 1-2j; 

Dj linear interpolation between vertices (Dj 
is D of Zi): and 

Lj ^ linear interpolation between vertices (Lj 
is L of Zj): 

di 2 x(Di*i +Di)xAZi; 
L ^ 2 x(U.i +Li)xAZi; 

where di and li are approximations to the L and D contributions for the i'^ cell and are computed by 
averaging the opacity and luminosity, respectively, that a face supplies, with the opacity and luminosity, 
respectively, that the last processed face supplies. Also. 

Adj (Dj + i-Di)xAzi; and 
Ali- (Li.i-Li)xAZi. 

The total transparency of a cell is given by: 



1 - dj ; or 

1 -di + ^xdiXdi; 



where the latter approximation is derived from a truncated Taylor series. 
The total brightness, taking into account opacity, is given by: 

li; or 
lixd-lxdi); or 

li X (1 - Ix di) + di X Adi - di X Ali); 

where the latter approximation is also derived from a truncated Taylor series. 

As was stated, the computation of the brightness may proceed from front to back, in which case for 
each segment and are updated by 
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T <- T Xtj 

5 

or the computation may proceed from back to front, in which case only X is needed: 

10 

While the front-to-back computation may appear to require more processing time, it has the advantage that 

the processing can be stopped when the accumulated transparency reaches a small enough value that 

further contributions are inconsequential. 
15 As indicated in Figs. 6 and 7 in relation to the regular grid compositing method of the prior art, a 

volume rendering computation requires accumulating, along the viewing ray that corresponds to each pixel. 

the light scattered from each piece of the volume is attenuated by passage through the volume between the 

light emitting portion and the viewer. Regular grid volume rendering is easily accomplished by 

"compositing" successive layers of the volume. 
20 By example, let U be the amount of light arriving at pixel r along the viewing ray associated with that 

pixel, and let l^j and Or.i represent, respectively, the luminosity and the opacity of the voxel where ray r 

intersects slice i. 

The well-known compositing algorithm for regular-grid volume rendering can be summarized as 
computing all the U by: 

25 

for each viewing ray r 

for each slice i of volume 
for each voxel v in the slice 
30 compute viewing ray r that intersects v 
compute lr,( and Or.i 

U-UX(I-Om) + lr.i 

However, as illustrated in Fig. 8 irregular grid volumetric data is typically represented as a three- 
35 dimensional mesh of space-filling polyhedra, or cells. Data values are known at the vertices of each cell, 
and can be interpolated from these vertex values to give data values throughout the volume of the cell. In 
the ray-tracing approach, volume rendering is accomplished by stepping along each ray from cell to cell, 
accumulating the light emitted by the volume as the ray passes through the cell, and attenuating the 
already accumulated light. 

40 By example, let L, be the amount of light arriving at pixel r along the viewing ray associated with that 
pixel, and let Irj and Or.i represent respectively the lunriinosity and the opacity of the ith cell of the mesh that 
ray r intersects. 

The known ray-tracing algorithm for irregular-grid volume rendering considers each ray r in turn and 
computes U for that ray by: 

45 

For each viewing ray r 

sort cells that ray r intersects 
for each cell i that ray r intersects , 
50 compute Or.i and of cell i 
U- U X (l-0,i) + Vi. 

In contradistinction, the object-based approach that is a feature of the Invention proceeds by computing 
U by rearranging the above computation so that the outermost loops consider each object in turn (hence 
55 object-based) in accordance with: 

sort cell faces 

for each viewing ray r 

6 
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for each cell face 1 

for each viewing ray r that intersects 

face i 

5 compute Or.i and l^.i of the cell 
behind face i 
U-UX(1-o,.i) + 

10 One significant advantage of the object-based technique over the ray-tracing technique is that the 
computation of the intersection between the rays and the cell faces may be rapidly accomplished. This is 
because each cell face is examined in turn and all ray intersections with the face are computed together, 
allowing for some portion of the intersection computation, such as a set up portion, to be shared between 
rays intersecting the same cell. A second advantage is that the general object-based approach is similar in 

75 some respects to the approach used in modern high-performance graphics systems for drawing conven- 
tional surface-based models, allowing this volume rendering approach to be more easily and completely 
integrated into such systems. For example, the ray-cell intersection computation may be accomplished in 
the same manner as that used for traditional surface rendering. The face-based orientation of the algorithm 
allows opaque and translucent surfaces, as well as volumes, all to be rendered together in the same image. 

20 Referring to Fig. 5 there is illustrated a face (F) of a tetrahedral cell having vertices Vi, V2, and V3. 
Associated with each vertex are a number of parameters including a depth (2) from the viewplane, RGB 
color information corresponding to the luminosity, and a value corresponding to the opacity (Op). Between 
vertices the data values are determined through a linear interpolation technique. 

Referring to Figs. 2 and 3 there is illustrated in greater detail the organization of the data structures, or 

25 lists LI. L2. and L3. In Fig. 2 it is seen that the list LI stores a plurality (1-n) of points corresponding to 
vertices of the faces of the grid cells. Each entry includes x, y. z positional data, relative to the viewplane. 
RGB color Information (luminosity) and the associated opacity (Op). L2 is a list of elements that specify the 
four vertices (V1-V4) for each tetrahedral cell (1-m). It should be realized that for other cell geometries more 
than four vertices may be required. Each of the four vertex entries for a given cell stores a vertex identifier 

30 and a pointer to one of the entries of LI. These pointers are indicated schematically by the arrows 
designated P in Fig. 2. As can be seen, some of the points in LI may be referenced by more than one 
vertex in L2, that is. certain of the vertices may have Identical positional, luminosity and opacity values. By 
this technique, the position, luminosity (RGB), and opacity of each cell vertex is defined and represented. It 
Is noted that in a conventional regular grid-based system, wherein all cell axes are orthogonally disposed 

35 relative to one another, that the list L2 Is not required in that the vertices are implicitly defined. 

Referring to Fig. 3 there Is shown a portion of the list L3 that is constructed by the algorithm during the 
process of rendering. By example, for a tetrahedral grid cell having vertices identified by (5. 8, 13, 22) the 
faces are defined by vertex identifier triplets (5, 8, 13), (5, 13, 22) etc.. The faces are ordered in the list L3 
by depth (z) from the viewplane. For the front-to-back method of processing the faces, the faces are 

40 ordered In L3 from closest to to farthest from the viewplane. For the back-to-front method of processing the 
faces, the faces are ordered In L3 from farthest from to closest to the viewplane. 

During operation, the processor 14 accesses the faces in L3 and processes them sequentially In 
accordance with the method described above to treat all rays that Intersect the face. While processing the 
faces the processor 14 constnjcts a displayable image in the image buffer IB. the displayable image 

45 eventually being provided to the frame buffer 18 for display upon the display 20. 

Although described in the context of a color display system it should be realized that the teaching of 
the invention also relates to black and white display systems. Also, it is within the scope of the invention to 
associate a separate opacity with each color for simulating the effect of preferentially attenuated light. 
Furthermore, the system embodiment disclosed in Fig. 1 is but one suitable embodiment for practicing the 

50 method of the invention. 

Thus, while the invention has been particulariy shown and described with respect to a preferred 
embodiment thereof, it will be understood by those skilled In the art that changes in form and details may 
be made therein without departing from the scope and spirit of the invention. 

55 Claims 

1. A method to be carried out by a computer system for rendering data for display, the data representing 
a model space, comprising the steps of: 

7 
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associating individual points in the nnodel space mth individual vertices of an irregular grid cell 
enclosed by faces, each of the points having an associated lunninesity value and an associated opacity 
value, whereby said step of associating preferably includes a step of determining an x, y and z 
5 coordinate of the point relative to the viewing plane and whereby said step of associating preferably 

includes a step of expressing the luminosity in terms of a plurality of colors; 

ordering faces of the grid cells, whereby said step of ordering preferably includes a step of ordering the 
faces by depth from the viewing plane; and 

10 

processing individual faces for determining, for alt viewing rays that project from pixels on a viewing 
plane and that intersect the face being processed, an amount of light arriving at the pixel, the amount of 
light being a function of both the luminosity and the opacity of a grid cell that lies adjacent to the face 
and through which the viewing ray passes, whereby said step of processing preferably is accomplished 
J5 by processing the faces in order of increasing or decreasing depth from the viewing plane. 

2. Apparatus for rendering data for display, the data representing a model space, comprising: 

means for associating individual points in the model space with individual vertices of an irregular grid 
20 cell enclosed by faces, each of the points having an associated luminosity value and an associated 
opacity value, whereby said means for associating preferably includes means for detemnining an x. y 
and 2 coordinate of the point relative to the viewing plane, and whereby said means for associating 
preferably includes means for expressing the luminosity in terms of a plurality of colors; 

25 means for ordering faces of the grid cells, whereby said means for ordering preferably orders the faces 
by depth from the viewing plane; and 

means for processing individual faces for determining, for all viewing rays that project from pixels on a 
viewing plane and that intersect the face being processed, an amount of light arriving at the pixel, the 
30 amount of light being a function of both the luminosity and the opacity of a grid cell that lies adjacent to 
the face and through which the viewing ray passes, whereby said means for processing preferably 
processes the faces in order of increasing or decreasing depth from the viewing plane. 

3. Apparatus as set forth in claim 2 and further comprising means, coupled to the expressing means, for 
35 displaying the colors. 

4. A method to be carried out by a computer system for rendering data for display, the data representing 
a model space, comprising the steps of: 

40 associating individual points in the model space with individual vertices of an irregular grid cell 
enclosed by faces, each of the points having an associated luminosity value and an associated opacity 
value; 

ordering faces of the grid cells by depth from a viewing plane; and 

45 

processing, in order, individual faces for determining, for all viewing rays that project from pixels on the 
viewing plane and that intersect the face being processed, an amount of light arriving at the pixel, the 
amount of light (L) being a function of both the luminosity (I) and the opacity (o) of a grid cell that lies 
adjacent to the face and through which the viewing ray passes, wherein 

50 

L— L X (1 - o) + 1. 

5. A method as set forth in claim 4, wherein the step of processing includes a step of estimating a total 
transparency (t|) of a grid cell in accordance with the expression 

55 

li 1 -di+ixdjxdi; 
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where Sj is an average opacity. 

6. A method as set forth in claim 4 or 5, wherein the step of processing includes a step of estimating a 
total brightness B| of a grid cell in accordance with the expression 

5 

where di is an average opacity and where Ij is an average luminosity. 

70 7. A method as set forth in any one of the claims 4 to 6, wherein X is an accumulated light at a pixel, 
where t is an accumulated transparency of processed grid cells, wherein the faces are ordered from 
closest to furthest from the viewing plane, and wherein X is updated in accordance with 

X ^ X + r X B: ; 

75 

and where r is updated in accordance with 



20 ^ T X t; - 

where Bj is an estimate brightness of a grid cell and where ti is an estimated total transparency of the 
grid cell. 

25 

8. A method as set forth in any one of the claims 4 to 7. wherein X is an accumulated light at a pixel, 
wherein the faces are ordered from furthest to closest to the viewing plane, and wherein X is updated in 
accordance with 

30 

X <- ^ X I: + 5: 

where Bj is an estimate brightness of a grid cell and where t j is an estimated total transparency of the 
35 grid cell. 

9. A method as set forth in any one of the claims 4 to 8, wherein the step of associating includes the 
steps of: 

40 maintaining a first list having entries including a three-dimensional position, a luminosity, and an opacity 
of points corresponding to vertices of the grid cells; 

maintaining a second list having entries including an identifier for each vertex of each grid cell; and 

45 establishing a pointer between each entry in the second list and an entry of the first list for specifying 
the position, luminosity and opacity of each vertex. 



50 



55 
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FIG. 3 
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® Object-based irregular-grid volume rendering. 
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@ Method and apparatus for rendering data for 
display. The data represents a model space and the 
method includes the following steps. A first step 
associates individual points in the model space with 
individual vertices (V) of an irregular grid cell en- 
closed by faces. Each of the points has an asso- 
ciated luminosity value and an associated opacity 
value. A next step orders the faces of the grid cells. 
The faces are preferably ordered by depth referen- 
ced to a view plane. A further step processes in- 
dividual faces for determining, for all viewing rays 
that project from pixels on the view plane and that 
intersect the face being processed, an amount of 
light arriving at the pixel. The amount of light arriving 
at the pixel is a function of both the luminosity and 
the opacity of a grid cell that lies adjacent to the 
face and through which the viewing ray passes. 
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